Non-universal Efimov Atom-Dimer Resonances in a Three-Component Mixture of 6 Li 
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We observed an enhanced atom-dimer relaxation due to the existence of Efimov states in a three- 
component mixture of 6 Li atoms. We measured the magnetic-field dependence of the atom-dimer loss 
coefficient in the mixture of atoms in state |1) and dimers formed in states |2) and |3), and found two 
peaks corresponding to the degeneracy points of the 1 23) dimer energy level and energy levels of Efi- 
mov trimers. We found that the locations of these peaks disagree with universal theory predictions, in a 
way that cannot be explained by non-universal two-body properties. We constructed theoretical models 
that characterize the non-universal three-body physics of three-component 6 Li atoms in the low energy 
domain. 



Since the first experimental evidence of Efimov states [Q]] 
in an ultracold cesium atomic gas |2|. few-body physics in 
ultracold atoms has attracted growing interest. The obser- 
vation of three -body and atom-dimer loss peaks and dips 
confirmed very general properties of few-boson systems 
near unitarity such as the universal scaling laws associated 
with the existence of Efimov states 10-01 ■ Although these 
loss enhancements and recombination minima are quali- 
tatively explained by the universal theory, their positions 
are shifted from universal predictions. While some efforts 
have been made to explain such resonance shifts by taking 
into account finite -range corrections ns MlQil . it is still an 
open question how accurately those corrections reproduce 
the observed Efimov spectra. 

Recently, it has turned out that a three-component Fermi 
gas of 6 Li offers another intriguing system to investi- 
gate universal few -body physics. Ottenstein et q/. llTlll and 
Huckans et al. lfl2ll observed the enhancement of the three- 
body loss at 130 G and 500 G in the mixture of fermionic 
6 Li atoms in the three lowest-energy hyperfine states 
\F;m F ) = |l/2;l/2), 1 1/2; -1/2) and |3/2; -3/2), 
which we label as |1), |2) and |3) respectively. Another 
resonance was later observed at 895 G fl 1 3M . It was ar- 
gued that those three -body loss enhancements are due to 
the existence of Efimov states Ill4l4l8ll . as in the case of 
identical bosons. In such a system, however, because of 
their fermionic nature, only distinguishable particles can 
interact via s-wave scattering. Thus interactions for the re- 
spective combinations are described by three different two- 
body scattering lengths a 12 , a 2 3 and a 3 i, which diverge at 
834 G, 81 1 G, and 690 G respectively due to Feshbach res- 
onances (see Fig. [Ha)). For each of these resonances there 
is a weakly bound dimer state which we designate as 1 12), 
1 23) and |31), respectively (see Fig. QIb)). Because these 
Feshbach resonance regimes overlap, the three scattering 
lengths can be varied simultaneously with a magnetic field, 
making it possible to realize a strongly-interacting multi- 
component system and access a whole new quantum phase 
of matter 119J. For this purpose, an accurate understanding 
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FIG. 1. (a) Two-body s-wave scattering lengths in 6 Li for states 
[1), 1 2), and |3) |20]. (b) Binding energies E h of the dimer 
and trimer states in three-component mixture of 6 Li. Dashed 
curves represent the energies from universal theory. Solid curves 
for dimers are obtained from two-body coupled-channel calcu- 
lations. Solid curves for trimers are obtained from our non- 
universal models adjusted to the experimental data. The inset 
shows the dimer-trimer crossing near 685 G. Dotted lines corre- 
spond to non-universal models without adjustment of the three- 
body parameter. 



of the low-energy few-body physics is necessary. As in the 
case of identical bosons, this requires an accurate determi- 
nation of the Efimov spectrum. 

In this Letter, we report the observation of two reso- 
nantly enhanced atom-dimer loss peaks in an atom-dimer 
mixture of 6 Li corresponding to the degeneracy points be- 
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tween the binding energy of 1 23) dimers and the ground 
(n = 0) and excited (n = 1) Efimov trimer states (see Fig. 
[Tib)). These peaks have been predicted by Braaten et al. 
from theory based on universality 111 811 . We compare the 
observed peak positions with universal theory predictions, 
and find significant deviations which cannot be explained 
by two-body physics only. We then construct non-universal 
models to interpret them. 

In our experiment we used an all-optical method to pre- 
pare a degenerate two-component Fermi gas of 6 Li atoms 
in the two lowest hyperfine states of 1 1) and |2) as de- 
scribed in detail in Bill . To prepare a mixture of |1) 
atoms and 1 23) dimers in equal population, we started 
with an imbalanced mixture of atoms in state |1) and |2) 
whose population ratio is |1) : |2) = 2 : 1. Evapora- 
tive cooling was performed at the magnetic field of 300 G 
where the amplitude of the scattering length between |1) 
and 1 2) shows a local maximum. The total number of 
atoms before dimer creation was ~ 10 6 . To achieve a 
very low collision energy, we adiabatically transferred this 
mixture into a larger volume hybrid magnetic/optical trap 
with smaller oscillation frequencies. The trap frequencies 
for atoms uj a were experimentally measured and approxi- 
mately given by w A = 2vr x VHO 2 - 0.735 Hz, w A = 

2vr x V50 2 - 0.315 Hz and u A = 2ir x \/7.2 2 + 0.455 
Hz in x, y and z directions respectively, where B is the 
strength of magnetic field in Gauss. Since the magnetic 
moments of 1 1) atoms and 1 23) dimers are different, the 
trap frequencies for dimers uP are different from the ones 
for 1 1) atoms, especially below 650 G. The trap frequencies 
for dimers were calculated using the magnetic moment of 
1 23) dimers. Deviation of uj b from the co 
directions at 580 G were about +6, +27 
tively. 

To create the atom-dimer mixture, we used a multiple- 
stage adiabatic rapid passage (ARP) as shown in Fig. E^a). 
By applying ARP at 563 G, we first transferred the atoms in 
state 1 2) to state |3). Then we transferred the atoms in state 
|1) to state 1 2) by ARP. Thus we created an imbalanced 
mixture of 1 2) and |3) whose population ratio is |2) : |3) = 
2 : 1. After the preparation of the imbalanced |2) - |3) 
mixture, we quickly ramped the magnetic field up to 8 1 1 G, 
and then swept down to the field of interest (580 G-760 G) 
in 300 ms to adiabatically create 1 23) dimers. At this point, 
excess atoms in state |2) were still left in the trap and the 
last stage of ARP was applied to transfer the atoms in state 
1 2) to state |1) to obtain a |1) atom - 1 23) dimer mixture. 

To study the decay of an atom-dimer mixture, we mea- 
sured the remaining fraction of atoms in 1 1) after a holding 
time of 1 ~ 1000 ms. The magnetic-field dependence of 
the atom-dimer loss rate was measured by repeating this 
sequence with various magnetic field values. Therefore the 
last stage of ARP was applied at a different magnetic field, 
and the condition for perfect ARP was checked at each 
magnetic field. When we take the absorption images, we 
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FIG. 2. (a) Experimental procedure to prepare atom-dimer mix- 
tures by using a multiple-stage adiabatic rapid passage (see text). 
Because the dimers cannot be imaged by resonance light, we dis- 
sociate the dimers by applying a magnetic field pulse to check the 
population balance of the atoms and dimers. (b) and (c): Typi- 
cal time evolution of the number of atoms in state |1) at 660 G 
and 595 G. The blue dotted curve shows an analytic fitting func- 
tion assuming the pure two-body loss and the black solid curve is 
obtained by a numerical fitting that includes the effect of dimer- 
dimer loss. Gray dashed curves show time evolution of the num- 
bers of molecules calculated from the rate equations. 



pulsed on the magnetic field gradient and spatially sepa- 
rated each spin component. This allowed us to check that 
there was no excess atoms in states |2) and |3) during the 
loss measurements. Although the temperature of the atom- 
dimer mixture is very low (~100 nK), it is close to the 
Fermi temperature due to the very low density. There- 
fore, we assume that the density distributions of atoms and 
dimers are Gaussian. Then, the number of atoms N A (t) 
and that of dimers N D (t) evolve in time according to the 
following coupled rate equations, 

N A = -FN A - CP^N A , (1) 

N D = -rW D - CP^N D - a^N u , (2) 

where f3, a and = 10 s are the atom-dimer loss 
coefficient, the dimer-dimer loss coefficient jind the one 
body loss rate, respectively. Here, Va 



and Vr, 



V A (0 A /V2u u f are the 



effective atomic and molecular volumes, where cr A , 



^k B T/m(u A y J 2 and a%. 



■ A ,,:) 2 and a» ytZ = ^ h !; T 2,„U n , ( ..) 2 
are the atomic and dimer cloud widths. We denote the 
mass of 6 Li by m, and u) A and up are the geometric 
means of the trap frequencies for atoms and dimers. C = 
EL \/4/(2+ {ujf/ujfY) ~ 8/\/27 is a numerical con- 
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FIG. 3. Magnetic-field dependence of the atom-dimer loss co- 
efficient f3 (red circles) and dimer-dimer loss coefficient a (blue 
diamonds) in the mixture of atoms in state |1) and dimers of |23). 
The black dashed curve is the calculated /3 from universal theory. 
Red, blue, and green solid curves show the /3 values calculated 
from our non-universal models A, B and B ' respectively. 



stant which results from the difference of the density dis- 
tributions of the atoms and dimers due to unequal masses 

a. 

We determined the atom-dimer loss coefficient j3 by fit- 
ting our data with the solution of the rate equations using V 
and a which were experimentally determined in advance. 
We checked that the cloud size remains approximately con- 
stant; therefore we assume the constant volume and ignore 
heating effect in our analysis. We found that the decay time 
scale of atom-dimer loss above 650 G is far shorter than 
the time scales of dimer-dimer and one-body loss and the 
decay curve can be described using two-body loss fitting. 
However, below 650 G, the dimer-dimer loss and one -body 
loss need to be taken into account (see Fig. |2b) and (c)). 

Figure [3] shows the measurement of the atom-dimer loss 
coefficient j3 (red circles) and the dimer-dimer loss coeffi- 
cient a (blue diamonds) for the 1 23) dimers, as a function 
of magnetic field. The error bars for a show the system- 
atic error of our measurements arising from the uncertainty 
in the number of atoms and size measurements. The error 
bars for (3 arise from the statistical error of size mesure- 
ments and error propagation of a. The systematic error 
arising from the uncertainty in the number of atoms, which 
is estimated to be ±22% and may shift the data linearly, is 
not included in these error bars. We observe two maxima 
of the atom-dimer loss coefficient near 602 G and 685 G. 

The atom-dimer loss is expected to increase at the mag- 
netic field where the energy level of an Efimov state inter- 
sects with the energy level of the 1 23) dimer II 1 811 - Since 
there is no other mechanism that causes an enhanced atom- 



dimer loss, the two peaks observed in Fig. [3] should cor- 
respond to such Efimov resonances. Because the lower 
peak is expected to be associated with a "ground-state" Efi- 
mov trimer, our finding provides the first experimental ev- 
idence for this "ground-state" Efimov trimer in the three- 
component mixture of 6 Li. 
We now analyze the data. The usual universal theory 
for this kind of three -body systems relies only 
on the 3 scattering lengths to describe the two-body interac- 
tions and a short-range three-body parameter Ae 171 , where 
A fixes the phase of the three-body wave function at short 
distance and rj > phenomenologically models losses oc- 
curing when 3 atoms come close. These free parameters 
were previously determined to be A = 0.885 1 and 
rj = 0.016 from the zero-energy resonance at 895 G ly- 
ing right in the universal region of large scattering lengths. 
Figures. (Hb) and [3] show the universal predictions for 
the trimer energy and atom-dimer loss coefficient using 
these previously determined parameters. While they quali- 
tatively explain the presence of the two resonances, there is 
significant disagreement. This is not surprising because the 
resonances investigated here are significantly away from 
the universal region. Near the resonance at 602 G, the scat- 
tering lengths are on the order of the range of the inter- 
actions (the van der Waals length £ vc iw ~ 60 a ) and the 
binding energy of 1 23) dimers differs by 57 % from the uni- 
versal form —h 2 /mal 3 . Near the resonance at 685 G, the 
scattering lengths are much larger, yet the binding energy 
still differs from the universal form by 8 %. If we follow 
the predicted universal trimer energy, this deviation sug- 
gests that the resonance point should be shifted to 685 G, 
which is interestingly very close to the observed peak loca- 
tion. 

To check this point and make a consistent analysis of our 
data, the minimal requirement is to design a theory which 
accurately describes the non-universal two-body physics, 
in particular the dimer energies. The most straightforward 
way to achieve this is to solve the three-body problem with 
zero-range interactions parametrized by energy-dependent 
scattering lengths. This can be easily implemented using 
the Skorniakov - Ter-Martirosian coupled integral equa- 
tions M, S IH. The scattering length a(k) should 
accurately reproduce the asymptotic two-body physics at 

i-2.2 

positive and negative energies E = . In particu- 

fz2 2 

lar, the dimer binding energy E h = — — is given by 
l/a(in) = k . The universal limit is retrieved for an 
energy-independent scattering length \/a{k) = 1/a. A 
realistic analytical expression near a Feshbach resonance 
can be derived from a two-channel model, for instance with 
separable Gaussian interactions ifioh . It leads to 
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where the first term contains the resonance parameters a^g, 
B , A(5), and u{B) which are fitted to reproduce accu- 
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rately the zero-energy scattering length, the effective range, 
and the last dimer binding energy - all obtained from a 
two-body coupled-channel calculation for each hyperfine 
configuration. The high-energy term f{k) is found to be 
/a(&) = kErf(kb) which behaves as k at large k. While 
the high-energy (short-distance) behaviour is irrelevant for 
the low-energy two-body physics, it does change the three- 
body phase at short distance. It is known however that 
two-body physics only cannot determine that phase in gen- 
eral Il22ll . Although we could adjust the high-energy two- 
body form to effectively set the 3 -body phase llOfl . we rely 
for that purpose on the three -body parameter he 111 , which 
appears as an upper bound of the integral in the STM equa- 
tions. To check that our theory does not depend on the two- 
body high-energy form, we consider an alternative form 
/b(&) = 1+ b (b fc )2 which tends to a constant | at large k. 
In all cases, we choose a form which is independent of the 
low-energy parameters, and is therefore magnetic -field in- 
dependent. 

Using these two models, which we refer to A and B, we 
first adjust the three-body parameters A and 77 to fit the 
previously measured loss coefficient near the zero-energy 
resonance at 895 G. As in the usual Efrniov theory, the 
choice of A is not unique and for completeness we consider 
two values for model B, which we refer to B and B'. We 
found that all models yield exactly the same universal res- 
onant loss profile near 895 G. Using those parameters, we 
then investigate the atom-dimer properties. All models give 
essentially the same prediction for the second resonance 
near at 685 G, but different ones for the first resonance at 
602 G, probably reflecting its strongly non-universal char- 
acter. The predicted location of the second resonance is 
around 672 G (see inset of Fig. [JJ, which is not shifted 
from the universal prediction as we naively expected - we 
also checked that the same result is obtained with a sim- 
ple effective range model. We conclude that our experi- 
mental data show that the short-range physics parametriza- 
tion cannot be constant altogether. Therefore, if we keep 
our two-body parametrization, the three-body parameter A 
must depend on energy, and possibly magnetic field. This 
seems reasonable, since the two-body parameters a{k) are 
already required to be energy-dependent in order to de- 
scribe the non-universal two-body physics. 

We proceed to map out this energy dependence by fitting 
each resonance with our models - see Table|H The variation 
of A with magnetic field is smooth and almost linear, but 
it corresponds to a non linear dependence on energy. The 
variation of ij is not monotonic and varies over one order of 
magnitude. This is not totally unexpected since it describes 
underlying loss processes which are at present unknown. 
We then perform a quadratic fit for both A and In r\ to es- 
timate those parameters at any magnetic field, and plot the 
corresponding atom-dimer loss coefficients for each model, 
as well as the Efimov trimer energies - see Figs. [Tib) and 
[3j The results are essentially the same, indicating that our 



TABLE I. Fitted three-body parameters A and 77 for different 
magnetic fields (energies) corresponding to different resonant 
points. These values are obtained for b — 38.85 clq in Model 
A, and b = 43.93 ao in Model B and B', which best reproduce 
the effective ranges of the interactions. We assume that the values 
for A are the same at zero energy. 



Model 


602 G 


685 G 


810 G(E = 0) 


895 G (E = 0) 


A 


A=0.0405 


A=0.0460 


A=0.0484 


A=0.0484 


B 


A=0.0111 


A=0.0152 


A=0.0186 


A=0.0186 


B' 


A=0.141 


A=0.185 


A=0.244 


A=0.244 


A 


77=0.008 


77=0.08 




77=0.033 


B 


77=0.020 


77=0.20 




77=0.011 


B' 


77=0.035 


77=0.30 




77=0.016 



analysis is model-independent. It reproduces the experi- 
mental data for B > 610 G up to a factor of 2. The first 
resonance cannot be very well reproduced, most probably 
because it is located at an energy which lies at the margin 
of validity of our zero-range calculation. 

In summary, we have measured the atom-dimer loss co- 
efficient (3 in a mixture of 6 Li hyperfine state |1) and 6 Li 
dimers in state 1 23), and found two loss maxima near 602 G 
and 685 G. We attributed these peaks to the crossings be- 
tween the atom-dimer threshold and two (ground and ex- 
cited) Efimov states. We found significant deviations from 
the universal predictions, and showed that they cannot be 
explained simply by taking into account the non-universal 
two-body physics. Our work therefore provides evidence 
for the non-universal character of short-range three-body 
behaviour, which we quantified by variations of the three- 
body short-range parameters. Although these parameters 
are different for different two-body models, they lead to 
a model-independent interpretation of our data, which pre- 
dicts the non-universal properties of the two Efimov trimers 
in three-component 6 Li. Understanding these three-body 
short-range variations will be a challenging task in the fu- 
ture. 

During the preparation of this paper, we became aware 
of similar results reported by T. Lompe et al. 12511 . 
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